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We study resistor diode percolation at the transition from the non-percolating to the directed 
percolating phase. We derive a field theoretic Hamiltonian which describes not only geometric 
aspects of directed percolation clusters but also their electric transport properties. By employing 
. renormalization group methods we determine the average two-port resistance of critical clusters, 

which is governed by a resistance exponent <j>. We calculate <fi to two-loop order. 

PACS numbers: 64.60.Ak, 05.60.-k, 72.80.Ng 
• Keywords: Directed percolation, transport, renormalization group, field theory, critical exponents 

o 



-a 



INTRODUCTION 



Percolation |Lj is a leading geometric model for connectivity in irregular media. In abstract terms, percolation is 
the passage of an effect through a medium which is irregularly structured in the sense that the effect can propagate 
through some regions whereas it cannot pass other areas. Such a medium could be a network of resistors connecting 
nearest neighboring sites of a d-dimensional hyper-cubic lattice and the effect would be a current / injected at site 
x and extracted at site x' . Irregularity in the medium could be realized by randomly removing resistors with a 
i-^ \ probability 1 ~ p. 

Depending on the occupation probability p the resistors (bonds) are likely to either be isolated or to form clusters. 
Two sites belong to the same cluster if they are connected by a path of bonds and hence current can flow between 
them. At low p two infinitely separated terminal sites x and x' are not connected by such a path and the network 
behaves as an insulator. For large p, on the other hand, many paths between x and x' may exist and the network is a 
y—i \ conductor. Therefore, at some probability in between, a threshold p c must exist where for the first time current can 
^ ■ percolate from x to x' . Below p c we have an insulator. Above p c we have a conductor. The threshold probability is 
called the percolation threshold, or, since it separates two different phases, the critical probability. 

In this picture the effect can percolate through occupied bonds in all directions. The resulting clusters are typically 
isotropic in space. Therefore ordinary percolation is referred to as isotropic percolation (IP). The linear extend of 
the isotropic clusters can be characterized by the correlation length £ ~ \p — p c \ " , where v is the correlation length 
y—\ ' exponent of the IP universality class. 

Directed percolation (DP) Q is an anisotropic variant of percolation. The bonds function as diodes so that 
the current can percolate only along a given preferred direction. Typical DP clusters are anisotropic and they are 
characterized by two different correlation lengths: £|| (parallel to the preferred direction) and £j_ (perpendicular to 
it). As one approaches the critical probability, the two correlation lengths diverge with the exponents v\\ and v±_ of 
\ the DP universality class. 

DP is perhaps the simplest model resulting in branching self-affine objects. It has many potential applications, 
including fluid flow through porous media under gravity, hopping conductivity in a strong electric field gj, crack 
propagation B, and the propagation of surfaces at depinning transitions ||. Furthermore, it is related to epidemics 
with recovery|]6| and self-organized critical models Q. 
• i-H , If the preferred direction is interpreted as time, a DP cluster represents the history of a stochastic process. In this 
interpretation the DP universality class is the generic universality class for phase transitions between an active and 
5-H \ an inactive absorbing state. 

The transport properties of IP have been studied extensively S. For example, the average two-port resistance 
Mr of a random resistor network (RRN), i.e., a percolation model in which the bonds are assigned a finite nonzero 
conductance, is known to obey in the vicinity of the percolation threshold the power law |l^| 

M R (x,x')~\x-x'f F/ " (1.1) 

with 4> lF being referred to as the resistance exponent for IP. 

Compared to IP relatively little is known about transport in DP. In this paper we study the average resistance 
in DP by employing the powerful methods of renormalized field theory. We find that this quantity is governed in 
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the vicinity of the critical probability by a resistance exponent (j> analogous to cf> lp . For instance, measured in the 
preferred, time-like direction between two terminal sites x = (x± = 0,t) and x' = (x' ± = 0,i') it scales like 



M R {x,x')~(t-t'f v * . 
We obtain that <f> is given to second order in the departure e = 5 — d from 5 dimensions by 
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Our field theoretic approach thrives on our real world interpretation of Feynman diagrams which has proven to be 
successful in studying transport in IP via RRN @ If, 11 S H 111' 0- A real 

-world interpretation in the same 

spirit applies to DP also. We show that the Feynman diagrams utilized in this paper can be interpreted as being 
directed resistor networks. On one hand this remedies the complexity of the field theory and makes it more intuitive. 
On the other hand it simplifies the actual calculations. As shown below, the basic task in calculating the average 
resistance in DP is to calculate the resistance of Feynman diagrams. The remaining steps are completely analogous 
to those well known from previous field theoretic approaches jl^. 19 to usual, i.e., purely geometric DP. 
Brief account of our results was given in |2l| . The purpose of this paper i s to present our work in some detail. 



The 

plan of presentation is the follo wing : SectionTn] comprises the model. In Sec. II A background on rando m res istor diode 
networks is provided. Section II B discusses Ohm's and Kirchhoff's laws for these n etwor ks. Section II C introduces 
a generating function for Mr and explains the involved replica approach. In Sec. [I D we set u p a field theoretic 
Hamiltonian. Section III contains the core of our renormalizati on gro up analysis. In Sec. [II A| the diagrammatic 



elements of our perturbation calculation are presented. Section III B is devoted to the real- world interpretation of 
the Feynman diagrams constituted by these elements. In Sec. Ill C we set up a renormalization group equation. Its 
solution provides us with the scaling behavior of Mr. In Sec. HID we compare our result for to related results by 
other authors. Section IV, finally, contains concluding remarks. Technical details are relegated to Appendices |A| and 



II. THE MODEL 



A. The random resistor diode network 



Our approach is based on a model which captures both, IP and DP, namely the random resistor diode network 
(RDN) introduced by Redner [^2], |23|, ^4|. A RDN consists of a d-dimensional hyper-cubic lattice in which nearest- 
neighbor sites are connected by a resistor, a positive diode (conducting only in a preferred direction), a negative 
diode (conducting only opposite to the preferred direction) , or an insulator with respective probabilities p, p+ , p_ , 
and q = 1 — p — p+ — p- . In the three dimensional phase diagram (pictured as a tetrahedron spanned by the 
four probabilities) one finds a non-percolating and three percolating phases. The percolating phases are isotropic, 
positively directed, or negatively directed. Between the phases there are surfaces of continuous transitions. All four 
phases meet along a multicritical line, where < r := p+ = p_ < 1/2 and p = p c (r). On the entire multicritical line, 
i.e., independently of r, one finds the scaling properties of usual isotropic percolation (r = 0). For the crossover from 
IP to DP see, e.g., Ref. |§. 



B. Kirchhoff's laws 



To be specific we choose n = 1/Vffl(l, ■ ■ ■ , 1) as the preferred direction. We assume that the bonds k^j) between 
two nearest neighboring lattice sites i and j are directed so that kuj) ■ n > 0. We suppose that the directed bonds 
obey the non-linear Ohm's law 

where Vb = Vj — V{ is the voltage drop over the bond between sites j and i and 7& denotes the current flowing 
from j to i. In the following we drop the subscript whenever it is save. The bond conductances cr;, are random 
variables taking on the values a, <j6 (V), aO (— V), and with respective probabilities p, p + , p_, and q. a is a positive 
constant and 9 denotes the Heaviside function. Note that the diodes are idealized: under forward-bias voltage they 
behave as "ohmic" resistors whereas they are insulating under backward-bias voltage. To avoid confusion, we point 
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out that the round brackets in Eq. (2.1) contain the functional argument of the bond conductance. Because the bond 
conductance depends on the voltage drop over that bond only via a Heaviside function and sign (Vb) — sign (J^) we 
may write <Jb (Vb) = <Jb 

Assume that an external current I is injected at x and withdrawn at x'. It is understood that x and x' are connected. 
The power dissipated on the network is by definition 



P = I[V X -V X >] . 

Using Ohm's law it may be expressed entirely in terms of the voltages as 

P = R + (x, xT 1 Wx - V x ,} 2 = X>*> ( V k) V k = p ({v}) 



(2.2) 



(2.3) 



The sum is taken over all current carrying bonds (the backbone) between x and x' and {V} denotes the corresponding 
set of voltages. R+{x. x') stands for the macroscopic resistance when / is inserted at x and withdrawn at x' . Similarly 
one defines R-(x,x') as the macroscopic resistance when / is inserted at x' and withdrawn at x. The two quantities 
are related by R+(x,x') = R—(x',x). From the power one obtains Kirchhoff's first law 



(J) (J) 



as a consequence of the variation principle 

d 
dVi 



The summation in Eq. ([2.4| ) extends over the nearest neighbors of i and Ii is given by Ii — I (Si tX — 8i,x')- 
Alternatively to Eq. (p.3|), the power can be expressed in terms of the currents as 



P = R + (x, x') I 2 = J2"H (4) 4 = P ({I}) , 



(2.4) 



(2.5) 



(2.6) 



with {/} denoting the set of currents flowing through the individual bonds. It is understood that cr& 1 Ijj = 

whenever ct& (/&) = 0. Kirchhoff's second law, saying that the voltage drops along closed loops vanish, can be stated 
in terms of the variation principle 



d 



aid) 



p 



{""}•') 



= o, 



(2.7) 



i.e., there are no independent loop currents 1^ circulating around a complete set of independent closed loops. 



C. Generating function 

Our aim is to calculate the average resistance between two ports x and x' . In this section we set up a generating 
function for it. We demonstrate that this generating function indeed serves its purpose and explain how the average 
resistance can be extracted from it. 

The average resistance we are interested in is precisely defined by 

M R (x,x') = (x + {x,x')R + {x,x'))c/(x+{x^'))c. (2.8) 

(■■■)c denotes the average over all configurations C of the diluted lattice. x+( x i x ') is an indicator function that takes 
the value one if x and x' are positively connected, i.e., if / can percolate from x to x' , and zero otherwise. Note that 
(x+(x, x'))c = (X- (x' , x))c is nothing more than the usual DP correlation function. 

We follow an idea by Stephen |2(| and its generalization to networks of nonlinear resistors by Harris J27j and exploit 
correlation functions of 

^ x (x) = exp(iX-V x ^ , A^O, (2.9) 
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as generating functions of Mr. In writing Eq. we switched to D-fold replicated voltages Vi 



, • • • , V[ D ^ ) and imaginary currents A^ = iZj — ^ A 



, ■ • • , Xf' I • The correlation functions 



rep 



(2.10) 



are given by 



where 



3 a=l 



(2.11) 



D 



a=l v 6 



(2.12) 



and Z is the normalization 



j exp 



-2 P (W) 



(2.13) 



Note that we have introduced an additional power term ^ £\ Vi- This is necessary to give the integrals in Eqs. (2.11) 
and (2.13) a well defined meaning. Without this term the integrands depend only on voltage differences and the 
integrals are divergent. Physically the new term corresponds to grounding each lattice site by a capacitor of unit 
capacity. Th e orig inal situation may be restored by taking the limit of vanishing frequency, lo —> 0. 



From Eq. ( 2.11 ) the benefit of the replication procedure becomes evident. It provides us with an extra parameter D 
which we may tune to zero. In this replica limit the normalization denominator Z~ D goes to one and hence does not 
depend on the configurations C anymore. Then the only remaining dependence on C rests in the power P appearing 
in the exponential in Eq. (2.11). In th e replica limit, therefore, we just have to average this exponential instead of 
the entire right hand side of Eq. (2.11). This average then provides us with an effective power or Hamiltonian which 
serves as vantage point for all further calculations. The effective Hamiltonian will be discussed in S ec. | II D| . 



Now we come back to the role of Eq. ( 2.11 ) as a generating function. The integrations in Eqs. ( 2.11 ) and ( 2.15 ) 
are not Gaussian due to the 6 functions. Therefore we em ploy the saddle point method. The saddle point equation 
is nothing more than the variation princ iple stated in Eq. (2.5). Thus, the maximum of the integrand is determined 
by the solution of the circuit equations (|2.4) . Provided that the condition I 2 a holds, we obtain 



G 



(x,x',X\ 



= ( exp 



A 2 

— 2 R + i x > x ') 



(2.14) 



J / c 



up to an unimportant multiplicative constant which goes to one in the limit D — > 0. Upon expanding G around 
A 2 = 0, 



G 



(x,x',X\ = (x+(x,x')) 



c 



A 2 

1- —M R (x,x') + 



(2.15) 



we learn that the average resistance may be calculated via 

-i d 



M R (x,x') = (x+(x,x')) 



c 



d 



(-A 2 /2) 



G 



(x,x',\j 



A 2 =0 



(2.16) 



i.e., G is indeed a generating function for Mr. 
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At this point a comment on the nature of A is appropriate. We work near the limit when all the components of A 



are equal and continue to large imaginary values. Accordingly we set 27 

A (Q) = i\ + £ (Q) (2.17) 
with real Ao and £ < - Q \ Yl a =i £^ = 0- The saddle point approximation in Eq. (2.14) may be justified by demanding 

|A |»1. (2.18) 

On the other hand one has 



2 , Pi 



Thus, one can justify the expansion in Eq. ( 2.15| ) by invoking the conditions 



A 2 < D- 1 and e 2 < 1 



Note that the replica limit D — > allows for a simultaneous fulfillment of the conditions (2. IS) and (2.20) 



(2.19) 
(2.20) 



D. Field theoretic Hamiltonian 



In this section we device a field theo retic Hamiltonian for resistor diode percolation. Starting point will be the 
effective Hamiltonian announced in Sec. II C . First of all, however, we need to take care of some further regularization 
issues. 

Since infinite voltage drops between different clusters may occur, it is not guaranteed that Z stays finite, i.e., 
the limit lim^^o Z D is not well defined. This problem can be regularized by switching to voltage variables t? taking 
discrete values on a £>-dimensional torus which we refer to as the replica space. The voltages are discretized by setting 
t? = At?fc, where At? = "&m/M is the gap between successive voltages, *&m is a voltage cutoff, k is a D-dimensional 
integer, and M a positive integer. The components of k are restricted to — M < < M and periodic boundary 
conditions are realized by equating k^ = fc( Q )mod(2M). The continuum may be restored by taking §m — > oo and 
At? — * 0. By setting §m — $oM, M = m 2 , and, respectively, At? = t?o/w, the two limits can be taken simultaneously 
via rn — > oo. Note that the limit D — > has to be taken before any other limit. Since the voltages and A are 
conjugated variables, A is affected by the discretization as well: 



A = AM , AA At? = n/M, 



(2.21) 



where / is a Z?-dimensional integer taking the same values as k. This choice guarantees that the completeness and 
orthogonality relations 



(2M) 



L_£ ex p(iA-#) 



^A,0mod(2A/AA) 



and 



(2M) 



mod(2MAiJ) 



(2.22a) 



(2.22b) 



do hold. Equation ( |2.22|) p rovides us with a Fourier transform between the t?- and A-tori. 
Now we revisit Eq. (2.11). The replica limit provides us with the effective Hamiltonian 



H lep = - In ( exp 



\ p (W 



-In 



U dab -H 



Ob ) exp 



where the distribution function of the bond conductances is given by 

/ K) - V S [c?b - <r] + P+ S [at - <t8(V)] 
+ p-6 [crb - aO(-V)] +qS [a k ] ■ 



(2.23) 



(2.24) 



The 5-functions are interpreted so that / °° dx8 [x] = 1. Carrying out the average yields 



(2.25) 



where we have introduced 



-ftr(^) = In < q + p exp 



D 

P+ J] exp [--0 (V Q >) tf^ 2 



D 

[Jexp -|6>(-tf( Q ))^ Q 



a=l 
)2 



(2.26) 



Now recall our choice for A in Eq. ( [2.17 ). Since A and $ are related via Ohm's law Eq. (2.1), we have to choose i? 
consistently: 



i? (Q) =^o + C (Q) 

with real i? and ( (q) , Ef=i = 0. We impose the conditions 

|t?o| » 1 , « D- 1 , C 2 « 1 



(2.27) 



(2.28) 



(-i?o) + 9 (i?o) exp (-|# 



>(i?o)+#Mo)exp(--tf 2 



Under these conditions we can write 

A ^ i/j = In </ /'exp 

+ P 

By doing a little algebra we recast Eq. ( 2.29| ) as 

K(tf)=6 (tf ) K+ (d) + 9 (-<?„) if_ (tf 

Here we dropped a term 

«(i?o)ln[l-p-p + ]+fl(-i?o)ln[l-p-p- 
which is not depending on the bond conductances. Moreover, we defined 



K ± (3) 



= In 



i P + P± ( cr ^ 

1 exp tr 

1 - P - P± ^ 2 



(2.29) 



(2.30) 



(2.31) 



(2.32) 



Note that i^i ^i?^ are exponentially decreasing functions in replica space with a decay rate proportional to a 1 . Next 
we expand K I § J in terms of the = exp I iA • 



A^O 

x [0(A o )£+(a)+0(-A o )£- (X)] , 

where we have exploited that ($o) = 9 (Ao)- ^A^j stands for the Fourier transform of K± 

1 



(2.33) 



(2.34) 
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Because K± 
as 



are rotationally invariant functions in replica space, their Fourier transforms can be Taylor expanded 



= T± + W±X + 



(2.35) 



where t± and w± ~ cr 1 are expansion coefficients depending on p and p±. These expansion coefficients satisfy 
T±(p,p+,p-) = Tzp(p,p_,p + ) and w±(p,p+,p-) = w^(p,p-,p+). Since K± (jCj depends only on the square modulus 

into two parts, one being even and the other being odd under A 



of A it is easy to decompose K ($b) i 1 

k (4) = J2 ^ (0 *-x 0') { \ [k+ (a) + 



-A: 



A^O 
1 



[fi(Ao) 



-Ao)] 



if 4- I A ) - if_ I A 



(2.36) 



Now we insert Eq. ( 2.36| ) into Eq. (2.25). We also carry out a gradient expansion in position space which is justified 
because the interaction is short ranged not only in replica but also in position space. We find 



1 r 

2 



ffrep = -2J2J 



K+ A ) +-KL A 



if 

V-x (i) 



^A (* 

- if _ 1 a 



(2.37) 



We proceed with the usual coarse graining step and replace the ip? (i) by order parameter fields tp^ (x) which 

inherit the constraint A 7^ 0. We model the corresponding field theoretic Hamiltonian TL in the spirit of Landau as a 
mesoscopic free energy and introduce the Landau-Ginzburg- Wilson type functional 



H = J d d x|-^^_ x (x) [r-V 2 +w\ 2 

A#0 

+ (0(A o )-0(-A o ))vV]v>x(x) 

+ 1 



J2 *1>-X ( X ) ^-X ( X ) ^A+A' ( X ) 
A,A',A+A'^0 



A>0 



(2.38) 



where we have discarded all terms irrelevant in the renormalization group sense. The parameter r is the coarse grained 
relative of t + + r_. It specifies the "distance" from the critical surface under consideration, w ~ cr -1 is the coarse 
grained analog of w + + u;_. The vector v lies in the preferred direction, v = tin. r, w, and v depend on the three 
probabilities p, p +1 and p_. v is zero if p+ = p_. We will see as we go along that our Hamiltonian TL describes in 
the limit w — > the usual purely geometric DP. Indeed TL leads for w — > to exactly the same perturbation series as 
obtained in (T|,|T|,|20). 



III. RENORMALIZATION GROUP ANALYSIS 

In this section we utilize TL and calculate the generating function G(x, x', A) by employing field theory augmented 
by renormalization. For background on these methods we refer to Ref . . We perform a diagrammatic perturbation 
calculation up to two-loop order. In the remainder we drop the regularization term proportional to lu for simplicity. 
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A. Diagrammatic elements 

Instead of working with 7i directly, we recast Ti. into the form of a dynamic functional |3C| [H], [32| . This strategy is 
convenient because it simplifies the following calculations from the onset. Assuming u ^ we introduce new variables 
by setting 



x { \ = n • x = vpt , tp = \v\ 1/2 s, g = \ v \ 1/2 g. 



(3.1) 



On substituting Eq. ( |3.l| ) into Eq. ( 2.38 ) we obtain the dynamic functional 



J = d d± Xl _ dtl - s_a (x_l, t) \p (t - Vi + w\ 



A#0 



+ (0(Ao)-0(-Ao))^]* x (xj_,i) 



() 



A,A',A+A'^0 



(3.2) 



p (r - V 2 + + (6 (A ) - (-A )) 
:G(x ±) t,A) = <5(x ± )5(i). 



where dj_ = d — 1. In Eq. (3.2) we dropped a term containing a second derivative with respect to t because it is 
irrelevant co mpa red to the retained term containing d/dt. 

From Eq. (3.2) we gather the diagrammatic elements contributing to our renormalization group improved pertur- 
bation calculation. Dimensional analysis shows that the coupling constant g has the naive dimension (4 — dj_)/2, i.e., 

d = d±_ + 1 = 5 is the upper critical dimension. The Gaussian propagator G fxx, i, X^j is determined by the equation 

of motion 

~) 

dt J 

(3.3) 

Solving the equation of motion is straightforward. For the Fourier transform G ^p, t, \ j of G ^xi, t, A^ we obtain 

G (p, i, X) = G+ (p, t, X) + G_ (p, t, X) , (3.4) 
where p is the momentum conjugate to and 

G± (p, t,\j =6 (±t) 9 (±A ) exp ^tp (t + p 2 + wA 2 



(3.5) 



For the perturbation expansion it is sufficient to keep either G+ ^p, t, Xj or G_ ^p, t, Xj and hence we discard 
G- (p, t, X) . 



B. Resistance of Feynman diagrams 

We may rewrite G + ^p, t, )Cj as 

G + (p,t,X) = 6»(t)6»(A )exp -tp (t + p 2 + w\ 2 

- 8(t)exp[-tp(T + p 2 )]5 X5 , (3.6) 

i.e., the principal propagator G+ ^p, t, )Cj decomposes into a propagator carrying A's (conducting) and one not carrying 

A's (insulating). This allows for a schematic decomposition of the principal diagrams into sums of conducting diagrams 
consisting of conducting and insulating propagators. In F ig. U\ we list the result of the decomposition procedure. 

As in our previous work on transport in IP jll], [l^, |l3|7[l4|, [l5[] , the conducting diagrams may be viewed as 
resistor networks themselves with conducting propagators corresponding to conductors and insulating propagators 




+2 



+ 



/ B 





FIG. 1: Here we list the decomposition of the primary two leg diagrams (bold) into conducting diagrams composed of conducting 
(light) and insulating (dashed) propagators. The listing extends up to two-loop order. Note that the conducting diagrams 
inherit their combinatorial factor from their bold diagram. For example, the diagrams A and B introduced below have to be 
calculated with the same combinatorial factor, namely |. 



corresponding to open bonds. The "times" t appearing in the conducting propagators correspond to resistances and 
the replica variables iX to currents. The replica currents are conserved in each vertex and we may write for each edge 
i of a diagram, Xi = A^ ^A, {k}^ , where A is an external current and {k} denotes a complete set of independent loop 

currents. The A-dependent part of a diagram can be expressed in terms of its power P: 



exp 



-pw 



exp 



pwP (X,{«}) 



(3.7) 



The real- world interpretation suggests an alternative way of computing the Feynman diagrams. To evaluate the 
sums over independent loop currents, 



exp 

{k} 



pwP 



(*,{*}) 



(3.8) 



we employ the saddle point method under the conditions discussed a t th e end of Sec. [II Q . Note that the saddle point 
equation is nothing more than the variation principle stated in Eq. (2.7). Thus, solving the saddle point equ ations is 
equivalent to determining the total resistance R({ti}) of a diagram, and the saddle point evaluation of (3.8) yields 



exp 



-R ({U}) pwX 2 



(3.9) 



where we have omitted once more multiplicative factors which go to one for D — > 0. A completion of squares in the 
momenta renders the momentum integrations straightforward. Equally well we can use the saddle point method which 
is exact here since the momentum dependence is purely quadratic. After an expansion for small w\ 2 all diagrammatic 
contributions are of the form 



/(p 2 ,i,A 2 ) = I P (p 2 ,t)-I w (p 2 ,t)pwX 2 + 



f[dt t \l-R{{ti})pwX 2 + ■■■] D(p 2 ,t;{U}) 



(3.10) 
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where D (p 2 , t; is a typical integrand as known from the field theory of DP |ll| [T|, |o|. Concrete examples for 
calculations of this type are given in Appendix |a]. 



C. Renormalization and scaling 



We proceed with standard techniques of renormalized field theory [£9| . The ultraviolet divergences occurring in the 
diagrams can be absorbed by dimensional regularization. We employ the renormalization scheme 

s^s = Z 1/2 s, t^t = Z~ 1 Z t t, (3.11a) 
w — > w = Z~ x Z w w , p — > I > = Z~ x Z p p , (3.11b) 
= Z-^Z^Z^G-^u 1 ^^/ 2 , (3.11c) 

where e = 4 — d± and p is the usual inverse length scale. The factor G e = (47r)~ di / 2 r(l + e/2), with T denoting the 
Gamma function, is introduced for convenience. Z, Z T , Z p , and Z u are the usual DP Z-factors, which can be found to 
second order in e in the literature []l9| , [20| . Thus, Z^, remains to be determined. In a two-loop calculation we compute 
the A-dependent parts of the diagrams displayed in Fig. 1 in dimensional regularization. The corresponding results 
are listed in App. [b]. From the two-point vertex function, constituted by the inverse Gaussian propagator and the 
diagrams given in Fig. 1, we extract Z w by minimal subtraction which gives 

9 m 2 r 5 71 
Ye ~ 72 



Z w = 1 



3 u 
87 



64 e 



l ln U 



(3.12) 



There are two points which we would like to emphasize. First, non-primitive divergencies cancel each other which is 
an indication for the renormalizability of the perturbation expansion. Second, for w — > the decomposition in Fig. 1 
collapses to the usual DP diagrams with the correct fore- factors, i.e., we indeed retrieve the perturbation series for 
purely geometric DP. 

The unrenormalized theory has to be independent of the length scale p~ l introduced by renormalization. In 
particular, the unrenormalized connected N point correlation functions must be independent of p, i.e., 



dp 



= 



for all N. Eq. (3.13) translates via the Wilson functions defined by 

d 



7... (u) 

(3(u) 

K (u) 
(w (u) 

Cp(«) 



M-H-ln Z.. 
op 



= M 



du 
dp 
<91nT 



= M 



= M 



dp 
dhiw 



= (-£ + 7 + 2 7p -lu)u. 



7 P - 7r 



dp 
dhip 



dp 



1-lp: 



(3.13) 

(3.14a) 
(3.14b) 
(3.14c) 
(3.14d) 
(3.14e) 



where the bare quantities are kept fix while taking the derivatives, into the Gell-Mann-Low renormalization group 
equation 



d d d d 

p— — I- p— — h tk— + wQ u 



, d N 



dp du dr dw 
x Gn f I *±,pt, wA 2 | ;t,u, p\ = . 
The renormalization group equation is solved by the method of characteristics. At the infrared stable fixed point 



(3.15) 



u* , determined by (3 (u*) = 0, the solution reads 

G N (jx ± ,pi, wA 2 j ;r,u,pj 

= P' N/2 G N ({lx ± ,l c ?pt, ;C W A 2 } ; Kt, u*, lp) 



(3.16) 
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where 7* = j(u*), k* =k(u*), (* = Q p (u*) , and (* = (w («*)• 

To derive a scaling relation for the correlation functions, a dimensional analysis remains to be done. It yields 



G N Ux±,pt, wA 2 } ;t,u,is\ 



u, 1 



(3.17) 



Equation (3.16) in conjunction with Eq. (3.17) now gives 

Gat ^|xj_,pi, wA 2 } ;t,u,/iJ 

= i^+^ n / 2 g n (hx ± ,i z pt, r+'^wP}; r 1 

Here 

'4 



7 = - 

2+c; 
1 

2- K* 





' 25 


f 


288 


e 


/1- 






12 I 




- / 1 






16 [ 



hi 



67 59 



288 144 
107 17 

288 ~ 144 



hi 



mil 



(3.18) 

(3.19a) 
(3.19b) 
(3.19c) 



are the well known the critical exponents for DP which have been calculated previously to second order in e |l9[ po| 
In addition, we introduced in Eq. ( [3.18 ) the resistance exponent 



"x(2-C) = l+24^ 



151 _ 157 
288 ~ 144 



In - 



(3.20) 



Equation ( |3.18 ) implies the following scaling behavior of the two point correlation function G = G2 at criticality: 

G(|x_L-x^|,i-t',wA 2 ) 
= l d ± +r >G(l\x_ L -x' x \,l z (t-t'),l- 



(3.21) 



where we dropped several arguments for notational simplicity. In the following we set x'j_ = and t', again for the 
sake of simplicity. The choice I — |xjj _1 and a Taylor expansion of the right hand side of Eq. (3.21) lead to 



G Ix.UujA' 



IX^I 1 ^-"/! 



l + wX^xA^ f w . 



XL 



(3.22) 



with fx and f Wt i being scaling functions which vanish for vanishing argument. Equally well we can choose I = t x l z 
which then leads to 



G 



(|xJ,t, w X 2 ) = t ^-v)l* h (\^L 



x |i + ™aV/^ 2 (^) + 



(3.23) 



where v\\ = v±z . and f W) 2 are other scaling functions which tend to constants as their argument is taken to zero. 
Exploiting Eq. (fTl|) we deduce that 



Mr ~ 



(3.24) 



if measured along the preferred direction. For measurements in other directions it is appropriate to choose a length 
scale L and to express the longitudinal and the transverse coordinates in terms of L: \x±\ ~ L and x\\ ~ L z ~ T . 
With this choice the scaling function f Wi ± reduces to a constant and we obtain 



M R ~ L^^ ~ T 0/l/ " 



(3.25) 
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D. Comparison to numerical data 



Now we compare our result for <j) to numerical results available in the literature. We are not aware of any numerical 
results for 4> itself. However, several authors determined the conductivity exponent t governing the increase of the 
macroscopic conductivity £ ~ (p — p c ) . The conductivity exponent is related to the resistance exponent via the 
scaling relation |23| t — if) + (d — l)v± — v\\. Redner and Mueller |}3| determined t in two dimensions by Monte 
Carlo simulations: t(d = 2) = 0.6 ± 0.10. Arora et al. [|4| did analogue and numerical simulations leading to 
t(d = 2) = 0.73 ± 0.10. Another value for comparison is t(d = 2) ps 0.7, obtained by Redner p3| from a real space 
renormalization group calculation. 



Crudely evaluating the e-expansion in Eq. (3.20) for small spatial dimensions leads inevitably to poor quantitative 
predictions. Therefore it is appropriate to improve the e-expansion by incorporating rigorously known features. We 
carry out a rational approximation which takes into account that t(d = 1) = 0. Practically this is done by adding an 
appropriate third order term. We obtain the interpolation formula 

t fa (l - |) (2 + 0.2083 e + 0.0604 e 2 ) , (3.26) 

which leads to t(d = 2) ps 0.8. 



IV. CONCLUDING REMARKS 



The field theory presented in this paper was tailored to describe electric transport on DP clusters. However, it 
is also suitable to study purely geometric aspects of DP. Our Hamiltonian Tt provides alternative means to derive 
the usual critical exponents of the DP universality class. For w — > it leads to the same perturbation expansion as 
obtained in Refs. (l|,|2(|. 

Our real-world interpretation of Feynman diagrams, developed in the context of RRN, proves to be a powerful tool 
in studying transport in DP. Here we saw that it allowed for an elegant and intuitive way of determining the average 
resistance Mr. However, the real-world interpretation is not only fruitful in studying Mr. It also turns out to be 
beneficial in investigating multifractality in DP ]37| and in calculating fractal dimensions of DP clusters [^l], Q . 

The resistance exp onent <f> is found to be larger than the corresponding resistance exponent for the random resistor 
network (RRN) Jl^ ||f|. This is intuitively plausible since long tortuous paths that contribute to the macroscopic 
conductance in RRN are suppressed in DP. 

Our result for </> is for dimensions close to five the most accurate analytic estimate that we know of. In two 
dimensions our result shows reasonable agreement with the known numerical results. It is certainly desirable to have 
more and firmer numerical data for comparison with our analytic result, in particular in three dimensions. We hope 
that this paper triggers further simulations of transport in DP. 
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APPENDIX A: CALCULATION OF THE CONDUCTING DIAGRAMS 



This appendix outlines the calculation of the conducting Feynman diagrams listed in Fig. [l] in terms of examples. 
For briefness we focus on those parts of the diagrams proportional to w\ 2 . For the remaining parts the reader is 
referred to the literature on the field theory of DP jy^, . 

We start with diagram A which stands for 



A 



dt 



E 

P B 



exp 



{-pt 



2r + 2p 2 



,(«-A) 2 ]} , 



(Al) 



where J is an abbreviation for (2tt) d± J d d± p. As argued in Sec. 1IIB the summation over the loop currents k can 
be carried out by determining the total resistance of the diagram. Since A is composed of two resistors of resistance 
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* in parallel its total resistance is given by 

R(t) = t/2. 

Thus, we find by expanding for small w\ 2 [cf. Eq. fl3.10|) 1 that A takes the form 



A = -wX 



2 pV r °° 



dt p— exp (— 2ptr) j exp (— 2ptp 2 ^ 



(A2) 



(A3) 



Carrying out the momentum integration is straightforward. We obtain 

,2 



2 i poo 

16 (4tt)<W 2 J yy ' 



K2 P9 2 T(2-d ± /2) 2 _ d±/2 
16 (4tt)^/ 2 



(A4) 



Expansion for small e = 4 — d± than finally leads to the result stated in Eq. (Bl). Diagram B can be calculated by 
the same means. One simply has to substitute R(t) = * for the total resistance. 
Now we turn to diagram H which stands for 



H 



A 4 r °° 



2 

x exp 



f dhdhdh I 

Jo JpM - - 

{-p(*i+*a) 



2t + 2p 2 + wk 2 + w(k - A) 2 | 



x exp | - pt 3 3r + p 2 + k 2 + (p - k)' 

+ WTT 2 + w(ti — 7?) 2 + w(H — A) 2 j. 



(A5) 



By applying the usual rules for determining the total resistance of resistors added in parallel or in series we find that 
the total resistance of H is given by 



R (ii, <2i h) 



*1 + t 2 , + 7}t 2 + 2t\t 2 + §*1*3 + §^2*3 

2ti + 2* 2 + §* 3 



(A6) 



The momentum integrations are easily carried out by completing squares. After these steps we arrive at 

:P5 4 1 



H 



2 (Att)^J 



dt 1 dt 2 dt 3 exp [-t (2ii + 2t 2 + 3f 3 )] 



1 t\ + *| + ^t 2 + 2*1*2 + |*1*3 + 1*2*3 



( 2 ^ /2 [2* 1 + 2* 2 + f* 3 ] 



3 + 



(A7) 



The remaining integrations can be simplified by changing variables: *i — > i*y, *2 — > it(l — x — y) and *3 — ► Irta. This 



gives 



H = -„A^ 4 r ( 4 -^)^- 



12 (4tt)^ 



d±/2 f l r l- x 

dx 
[i Jo 



1 , 1 



y +-Ai-x-y) + 



~y(l - x - y) + ^a;y + ~z(l - a; - y) [> . 



(A8) 



where we have already carried out the integration over *. The integration over y can be looked up in a table ]3q |. We 
obtain 



H 



W r(4-d ± ) d| _ 4 /3 
~ WX "12" (4tt)^ T U 



dj_/2 



(A9) 
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where 



h = 

h = 
h = 



dxx 2 ~ d ^ 2 [1 + x] d± - 4 (1 - X) , 
dxx 1 -^/ 2 [1 + x] d± - A (1 - xf , 
dxx~ d±/2 [1 + x] d ^ i (1 - xf . 



(A10) 
(All) 
(A12) 



The integral I\ is convergent for d± < 4. Hence it is sufficient to evaluate it at the upper critical dimension, i.e., for 
dj_ = 4. Then we immediately obtain that 



h = 1/2 • 

Integral 1% can be evaluated by separating its divergent and convergent contributions: 



(A13) 



dx 



{ x l-d±/2 +x -l [ ( l_ a .)2_ 1 ]2| 



4-di 



I3 may be treated by the same means as I2 giving 



2-di 



2(d±-7) 5 
4 — d 1 2 



(A14) 



(A15) 



After substitu ting Eqs. (A13) to ( A15| ) into Eq. ( |A9| ) we carry out an expansion for small e which leads to the result 
stated in Eq. (B7). Diagrams I to L can be evaluated in an analogous fashion. 
As a third and final example we consider diagram C standing for 



C = P y 

x exp 
x exp 
x exp 



dtidt2dta 



k S 



{-fti 

{-ph 



2t + 2p 2 + wk 2 + w(k - A) 2 I 



2t + 2k + tun 2 + w(tt - 



A) 2 ]} 



3r + p 2 



(p-k)^ 



+ wit 2 + w(k — if) 2 + w(k — A) 2 I 
Solving Kirchhoff 's equations for this diagram results in 

1 

R (ti, *2, ^3) 



(A16) 



X |*i [2*!* 2 + 2*1*3 + 3*2*3 + 2*3] 2 

+ * 2 [2*1*2 + 2*2*3 + 3*1*3 + 2*2] 2 

+ (*i + * 3 ) [2*!* 2 + 2<i*3 + *2*3 + 4] 2 

+ (* 2 + * 3 ) [2*1*2 + 2*2*3 + *1*3 + 4] 2 



+ *3 [*l*3 + *2*3 + * 2 ] 2 



Hence the summations over the loop currents lead in conjunction with the momentum integrations to 

C = -w\ 2 pg 4 {M 1 + 2M 2 + 2M 3 } , 



(A17) 



(A18) 
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where 

Mi 



(47 



1 Z" 00 

—5- / dtidt 2 dh exp -t (2t a + 2i 2 + 3i 3 ) 
n) a± Jo 

[ti+t 2 + t 3 ] 2 

Wi + h){t 2 + h)-tl 
l r°° 

M 2 = / dtidt 2 dH exp -T(2ti+2ta + 3* 3 

[2^2 + 2jit 3 + 3t 2 f 3 +2t|] 2 
[4(t! +t 3 )(t 2 +i 3 ) -*§] 



Ma 



(4 



1 [°° 

—j- / dt 1 dt 2 dt 3 exp -t (2ti + 2* 2 + 3t 3 
<T X Jo 



, 2 



2*it 2 + Mitg + fgtg + 4\ 
x (ti+ts) 1 —; -i-- (A21) 

[4(tl+t 3 )(*2+*3)-^] 

The integral Mi can be simplified by the following change of variables: t% — > i(a; — 1), t 2 — * i(y — 1), and £3 — > i. 
Upon carrying out the integration over t we then obtain 

Mi = T 4^r^ r dxd y 



[2x + 2y-l] d -- 4 [z + y-1] 



2 



[Axy - 1] 



2+dj_/2 



r(4 - d A 



_dj--4 -7 



dx / dy 



[«-$r* ■ ( ' 

Since the remaining integrations are convergent for d± < 4 we evaluate them directly at d± = 4 which gives in 
e-expansion 

M 1 = -^r-« {25 - 1081ng)}. (A23) 

M 2 and M3 can be analyzed similarly. However, these calculations are somewhat more tedious. After all one obtains 
for C the result stated in Eq. (B3). Diagrams D to G may be computed like C. 
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APPENDIX B: RESULTS FOR THE CONDUCTING DIAGRAMS 



Here we list our results for the conducting diagrams. For briefness we display only those parts of the diagrams 
proportional to wX 2 : 



A 
B 
C 
D 
E 

F = G 
H 
I 
J 



\*2 2 — e/2 

-WA PQ — — T 1 



G, 



- w X 2 pg 2 ri r" e / 2 
4e 



\ 2 4 

-wA pa — r 
o e 

T2 4 u £ 

-wX pg —t 
16 e 

r* 2 

-wA pg —t 
o e 

T/2 4 - 

Hy 8e 



G 



2 V 3 



= wA P5 



o ^2 
2 „„4 ° ^-e 



128 £ 



G 2 

wA 2 ^/ — - t" 
16 e 



w\ 2 pg 4 



128 e 

K = L = wAV^r- 
64 e 



1 1 1 /4 
7 + 2-2 ln U 
2 

- + 1 - 111 . ., 
e V 3 



1 55 1, /4 

— I h - In - 

e 36 2 I 3 



1 1 /4 

7 + 1 + 2 ln U 
7 + 5 + ln 3 



2 , /4 

7 + 5 + ln 3 



(Bl) 
(B2) 
(B3) 
(B4) 
(B5) 
(B6) 
(B7) 
(B8) 
(B9) 
(BIO) 
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